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ABSTRACT 

We calculate the heating and cooling of the accreting white dwarf (WD) in the ultracompact AM Canum 
Venaticorum (AM CVn) binaries and show that the WD can contribute significantly to their optical and ultra- 
violet emission. We estimate the WD's effective temperature, T^^, using the optical continuum for a number 
of observed binaries, and show that it agrees well with our theoretical calculations. Driven by gravitational 
radiation losses, the time averaged accretion rate, (M), decreases monotonically with increasing Poib, covering 
six orders of magnitude. If the short period (Porb < 10 min) systems accrete at a rate consistent with gravita- 
tional radiation via direct impact, we predict their unpulsed optical/UV light to be that of the Tgff > 50,000 K 
accreting WD. At longer Porb we calculate the Teff and absolute visual magnitude. My, that the accreting WD 
will have during low accretion states, and find that the WD naturally crosses the pulsational instability strip. 
Discovery and study of pulsations could allow for the measurement of the accumulated helium mass on the 
accreting WD, as well as its rotation rate. Accretion heats the WD core, but for P^^ > 40 minutes, the WD's 
Teff is set by its cooling as (M) plummets. For the two long period AM CVn binaries with measured parallaxes, 
GP Com and CE 315, we show that the optical broadband colors and intensity are that expected from a pure 
helium atmosphere WD. This confirms that the WD brightness sets the minimum light in wide AM CVn bi- 
naries, allowing for meaningful constraints on their population density from deep optical searches, both in the 
field and in Globular Clusters. 

Subject headings: binaries: close — gravitational waves-novae, cataclysmic — variables — white dwarfs 



1. INTRODUCTION 

The AM Canum Venaticorum (AM CVn) class of ultra- 
compact binaries provides an opportunity for probing the 
accretion-induced heating of a white dwarf (WD) over six or- 
ders of magnitude in the time-averaged accretion rate, (M). In 
these very tight (Porb < 80 min) binaries (see Warner 1995 for 
an overview), the donor star is a very low mass (M < O.IMq) 
helium WD and the accretor is either a He or a C/O WD. 
Though there are a few distinct AM CVn formation scenar- 
ios (e.g. Nelemans et al. 2001 and references therein), they 
all evolve as mass-transferring binaries driven by the loss of 
angular momentum from gravitational radiation (e.g. Deloye, 
Bildsten & Nelemans 2005). 

This paper is an outgrowth of Townsley and Bildsten's 
(2004, hereafter TB) work on the reheating of WDs in H- 
accreting cataclysmic variables. However, the case of WD 
reheating in AM CVn binaries is different because the ac- 
creted material is pure helium, making flashes an isolated 
rather than constant occurence, and the binary evolution pro- 
ceeds only under angular momentum loss by gravitational ra- 
diation, leading to a plummeting (M) with increasing PqA- 
We have found that the WD is heated in the earliest phases 



of binary evolution and cools at later times. There is a criti- 
cal decoupling accretion rate of Mjec ~ (1-3) x 1O~"*M0 yr"' 
below which the reheated WD will cool even while accreting. 
We find the amount of time that has passed since the binary 
had this (M) by using the evolutionary calculations of Deloye 
et al. (2005). Simple cooling after this point predicts the WD 
luminosity, as well as the expected WD Teff and absolute vi- 
sual magnitude. My (assuming, of course, a WD mass, M, 
and radius, R). We show that the accreting WD can become 
the dominant broadband optical light source in those binaries 
with orbital periods in excess of « 40 minutes. 

In §2, we discuss the physics of the WD heating due 
to accretion and discuss the uncertainty in the total helium 
layer mass expected on the WD. Since the WD luminos- 
ity, L, is nearly always less than the accretion luminosity, 
L-dcc = GM{M)/R, we explain in §3 how the binary's behav- 
ior allows the WD to be directly detected at different orbital 
periods. At early times, during direct impact accretion (e.g. 
Marsh and Steeghs 2002), the unpulsed UV/optical light is 
from the hot WD, whereas at longer orbital periods, the disk 
is often quiescent, allowing for the underlying WD to domi- 
nate the optical continuum of a number of these binaries. We 
also speculate on what could be learned if one of these WDs 
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was caught in the DB pulsational instabihty strip. We make 
direct comparisons to the observations in §4, and conclude in 
§5. 

2. THE PHYSICS OF WHITE DWARF HEATING IN AM CVN 
BINARIES 

The ability for accretion to heat the deep interior of a WD 

has been discussed extensively (Sion 1991; TB). Indeed, ac- 
cretion of pure helium at rapid rates M > 10"^ Mq yr"' is 
known to flash unstably and was studied as a Type la progen- 
itor scenario (e.g. Nomoto & Sugimoto 1977; Nomoto 1982; 
Woosley & Weaver 1994). Though eventually ruled out as a 
major mechanism for Type la supernovae (Solheim & Yungel- 
son 2005), AM CVn binaries accrete pure helium for a long 
enough time to reheat the accreting WD. An additional feature 
of these binaries is the large dynamic range in the time aver- 
aged accretion rate, (M), onto the WD, allowing us to test the 
theories of re-heated WDs, as discussed in §3 & 4. 

Here we calculate the accre tor's thermal evolution along 
with that of the binary for several AM CVn systems that rep- 
resent the expected ranges for this population. Motivated by 
the expected mass distribution of He WD donors that evolve 
into AM CVn systems (Nelemans et al. 2001), we chose an 
initial donor mass of O.24M0. The thermal state of the donor 
is very important for setting the relationship between (M) and 
Porb (Bildsten 2002; Deloye & Bildsten 2003; Deloye et al. 
2005). For this study we chose two thermal states found by 
setting the initial entropy of the donor so that the system's Poib 
at contact is greater than the forb at contact of 30% and 90% 
of the systems in the Deloye et al. (2005) AM CVn popula- 
tion evolution calculation. This provides (M) as a function 
of time, and we then use the quasi-static envelope methods of 
TB to heat the accreting WD. For simplicity, this evolution 
is performed at constant M and with a constant helium layer. 
We chose two masses M = 0.65Mq and 1.05Mq, each case 
with a O.O5M0 accreted heUum shell (justified below). Initial 
M values for the binary evolution were chosen such that M(t) 
matched that of the thermal calculations at late times (long 
Porb) to ensure the correct Porb-time relation. 

Using these (M) trajectories (see Figure 1), we identi- 
fied four relevant regimes for the state of the accreting WD. 
Chronologically since initiation of mass transfer (and in or- 
der of decreasing (M)), these are: (1) rapid accretion and 
thermally unstable He ignition (< few x lO^yr), (2) core heat- 
ing during a phase in which L is determined from "compres- 
sional" heating in the WD envelope ('--^ 3 x 10^ - 3 x lO^yr), 
(3) a transition period where the core temperature is near an 
equilibrium value set by (M) (3 x 10* - lO^yr), and (4) the 
WD cooling regime (> 10^ yr) where L is determined by the 
time, Af , since the heating stopped at a previous epoch where 
(M) ^ Mdec, the decoupUng mass-transfer rate. We now de- 
scribe each regime. 

At the very earUest times, when (M) > IQ'^Mq yr"\ the 
accreted helium is compressed rapidly and ignites in a ther- 
mally unstable manner, leading to recurrent weak explosions 
(Fujimoto & Sugimoto 1982).' It is uncertain how much of 
the accreted material is ejected (see Kato & Hachisu 1999). 
As (M) drops, the compression happens more slowly so that 

' Extremely liigli accretion rates for which a red giant envelope can 
be supported, (M) > Lhc Edd /cn = 7 X Kr^Mg yr"'(M/MQ) (where c„ = 
0.6 MeV/amu is the energy released from helium fusion) only appear at the 
brief (< 10' yr) and poorly understood onset of accretion. 




Fig. 1 . — The time averaged mass transfer rate (M), WD core temperature, 
Tc, and WD luminosity, L, as a function of time since the onset of mass trans- 
fer for donors with fixed entropy. The thin (thick) lines show the evolution for 
a 0.65Mq (I.OSMq) WD, 0.05Mq of which is a surface helium layer. For 
the lower mass we show two curves with low (solid, cool) and high (dashed, 
hot) donor entropy (from Deloye et al. 2005). There are four epochs: (1) (not 
shown) rapid accretion with heUiun shell flashes (2) Tc rises at early times 
when (M) becomes low enough that shell flashes are suppressed and the ac- 
creted helium layer can build up, (3) (M) falls enough that T^eq Ri Tc and 
core heating ceases, (4) (M) becomes low enough that the cooling of the core 
dominates the outgoing luminosity. 



the envelope temperature is lower, and the heUum shell can 
become quite large. Thermal instabilities at this stage likely 
lead to detonations (Fujimoto & Sugimoto 1982), such that 
the last thermonuclear instability might be the most violent. 
However, once (M) < 2 x 10~'^ Mq yr"', the required unsta- 
ble He ignition mass for O.6-I.IM0 WDs (Nomoto 1982; 
Limongi & Tomambe 1991; Goriely et al. 2002) exceeds the 
donor mass. Thus once a system reaches an orbital period 
of 10 minutes, it is unlikely to undergo further unstable he- 
lium ignition events. That is, (M) drops quickly enough that 
once (M) is low enough for a strong flash, the accreted layer 
never becomes thick enough for it to occur. Over the remain- 
ing course of the binary's evolution, about O.IMq is accreted 
by the WD, motivating us to use a fixed accreted layer mass 
of O.O5M0 to calculate the thermal evolution of the accretor. 
Because L depends only logarithmically on the accreted layer 
mass (TB), using this fixed value has only a small impact, es- 
pecially when compared to the large variations in L from the 
orders of magnitude variation in (M). 

The period of time during which (M) > 1O"^M0 yr~' is so 
short (^ 10^ yr) that the accretor's core is thermally unper- 
turbed from its initial state. However, once the hehum layer 
begins to accumulate, thermal diffusion both out of the sur- 
face (setting the exiting luminosity, L) and into the core (act- 
ing to increase the temperature there, Tc) occurs. We thus 
begin evolving the WD's thermal state 10^ yr after contact. 
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Although the compressional energy release depends weakly 
on the accumulated He layer mass (i.e. logarithmically, TB), 
the L-Tc relation which determines the cooling is essentially 
independent of it. The evolution is shown in Figure 1 and ex- 
hibits all of the epochs (2)-(4) including reheating, decoupUng 
and cooling. The heat capacity and crystallized fraction of the 
core was calculated using the equation of state of Potekhin & 
Chabrier (2000), and the late-time L-Tc relation (in the limit 
of minimal accretion) was adapted from Hansen (1999). We 
have included a latent heat of crystallization of Q.llkT per 
particle, but no fractionation. 

At early times in the evolulion, Ihe WD core temperature, 
Tc, is much lower than the maximum temperature reached at 
the base of the outer radiative envelope, T^. In this limit, 
TB showed that L becomes independent of and is only 
a function of (M) and M. They found that the envelope 
state is determined self consistently by balancing the heat 
released in the outer (non-degenerate) envelope due to fluid 
motion down the entropy gradient, Lcomp ~ 2{M)kTi„/ fimp, 
with that carried through the envelope by radiative transfer 
L PS L©(r„/10'* K)^^, where /x = 4/3 is the mean molecular 
weight in the He envelope, k is Boltzmann's constant, and 
ntp is the baryon mass. This same state applies here when 
Tc < Tm and we have calculated this balance numerically us- 
ing the quasi-static envelope methods outlined in TB for He 
accretion and found. 
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This relation holds for 
10^ K((M)/lO-''M0yr-i)''-5 and IO-'^Mq yr"' < (M) < 
1O~^M0 yr"', although for our time histories Tc ~ T^ by the 
time (M) ~ 10~^ Mq yr"' so that this form only apphes to 
Figure 1 for (M) > lO'^M© yr"'. During this early phase, the 
WD core is being heated as (M) is dropping.-^ 

As (M) drops and Tc increases, the accreting WD ap- 
proaches the equilibrium state where the outgoing luminosity 
is matched by the internal energy release (Nomoto 1982; TB). 
There is a brief transition from core heating to cooling which 
appears in Figure 1 as a plateau in Tc. In the equilibrium state, 
a fairly good estimate of the compressional heating luminosity 
is found by integrating throughout the O.OSMq accreted de- 
generate He layer (identical to that performed in Appendix B 
of TB, and Nomoto 1982), giving Lcomp ~ lO{M)kTc/iJ.iflemp, 
where /i, He = 4 is the ion mean molecular weight in the he- 
hum envelope. As (M) continues to drop, this will eventually 
become low enough that it is exceeded by the core cooling lu- 
minosity Lcooi = Lq (Tc/ 10^ Kf -^ . At this point the WD evolu- 
tion decouples from (M) and cools much as an isolated object 
would. By setting Lcomp = Lcooi and taking Tc = 1-2x10^ K, 
we find Mdec = 1-3 x 10~^^Mq yr"'. These values are evident 
in Figure 1 . 

The precise value of Tc max at the plateau depends on the 
initial Tc and when the He flashes cease, neither of which is 

^ Tlie only otlier calculation we could find that reported exiting luminosi- 
ties are in Table 2 of Limongi & Tornambe (1991). While our value com- 
pares well to their calculations with M = 0.6Mq as the starting WD mass, 
when M = ().8Mq and the accreted layer is thin, their values are as much as 
a factor of 50 lower for (M) ~ Ur^A/Q yr"' . We believe that L, which is 
largely due to compressional heat release in the outer layers and was only of 
peripheral interest in their study, is not well determined by their modeling for 
small accreted layer masses at M = 0.85Mq . 



well known. This leads to some uncertainty in L at the time of 
decoupling. However, analogous to isolated cooling WDs, the 
"initial" Jcmax is forgotten, and L eventually only depends on 
the time after decoupling, At. This is simply the WD cooUng 
age from the time that heating ceased. By lO' yr, L » 1O~^L0, 
already consistent with a cooling DB WD of a similar age 
(Hansen 1999). 

Using the values tabulated by Limongi & Tornambe (1991, 
hereafter LT91) for accretion of Helium onto a O.8M0 WD 
at a fixed rate, we find that our determination of the heating 
time from the quasi-static model agrees well with their full 
numerical model. Table 2 in LT91 gives a detail of Tlnax (the 
maximum with respect to radius) and Tcenter for a mass history 
starting from O.8M0 up to about O.9M0. We compare our 
Tcore = Tc(t) for a O.9M0 WD with a O.O5M0 Helium layer 
accreting at a constant rate to T^ore = (7Lenter + 7max)/2 (repre- 
sentative of the heat content) from Table 2 of LT91, and find 
that in both cases, for this range of mass, dTcom/dM is approx- 
mately constant. We find dTcort/dM w 1.6 x 10^ K Mg' and 
2.7 X 10^ KMgi for (M) = 1O-^M0 yr"' and 5 x 1O-*'M0 yr"' 
respectively, compared to the values of 1.3 x 10** K Mq and 
2.7 X 10^ K from LT91. Also, the T^max reached in the 
plateau (^ 4 x 10^ K) reflects the equilibrium Tc for the (M) 
at this epoch ~ 5 x 1O~^M0 yr~'. This is consistent with the 
typical values of interior temperatures, log T ~ 7.5-7.7, found 
by LT91 for a O.8M0 CO WD accreting HeUum steadily at 
(M) ~ 5-10 X 1O-'^M0 yr-^ 

Before core crystallization begins, the WD's heat capacity 
is well approximated by C = 3kM / fiinip. With the simple re- 
lation L w Lq(Tc/\Q^ Kf \M/Q.6MQ), L = -CdTc/dt gives 

L w 6 X \G~^LQtg^^^{M /0.6Mq) where ^9 is time in Gyrs. This 
agrees well between decoupling and the onset of crystalliza- 
tion, which occurs at Tc = 4(9) x 10*' K forM = 0.65(1 .05)Mq, 
corresponding to a time of 1 .4(0.4) x lO** yr and an orbital pe- 
riod of 60(45) minutes for the hot donor evolution. The effect 
of the latent heat release is visible in the Tc history and the 
heavier WD cools more quickly at late times due to its earlier 
crystallization and the resulting decrease in total heat capac- 
ity. The modifications due to latent heat release change the 
L-t relation to something closer to L = 6.8 x lO'^^Lgfg ' ' and 
1.4 X IQ-^L©?,!-^ for M = 0.65Mq and l.OSM© respectively 
in the 1-8 Gyr region. By contrast, the accretion luminosity 
Lace = GM{M)/R is a smooth power law after about 10^ yr, 
given by L^c = 8.3(19) x IQ-^LQtg^-^ forM = 0.65(1. 05)Mq. 

The bottom panel of Figure 1 shows the resulting L(t), for 
which there is very little contrast with the donor entropies. 
However, when displayed as a function of the observable, PqA 
(see Figure 2), a contrast immediately appears: accretors in 
systems with hotter donors have higher Tiff's. The reason 
for this depends on the phase of evolution, but is essentially 
related to the difference in speed of evolution as set by the 
donor's M-R relation. In the "compressional heating" phase a 
lower entropy donor always produces a lower (M) at a given 
Porb, leading to a lower L. By an orbital period of 30-40 min- 
utes the WD is cooling and the time since decoupling becomes 
larger than the uncertainty introduced by the unknown Tc^ax, 
so that at longer orbital periods we can safely predict L. Low- 
entropy donors reach the decoupling (M) at shorter orbital 
periods and evolve more slowly in Porb, so that the WDs at a 
given Porb always look "older" (i.e. colder) if their donor was 
initially of low entropy. The track for the cool donor (30th 



4 



Bildsten, Townsley, Deloye and Nelemans 



M=1.05Mg + hot 
M=0.65Mg + hot -| 
M=0.65 M„ + cool 




10 20 30 40 50 60 70 80 90 

P^^i^ (minutes) 

Fig. 2. — The time-averaged accretion rate, (M), WD surface tempera- 
ture, Teff. and minimum absolute visual magnitude. My for the accreting WD 
in AM CVn binaries. The different WD curves are as in Figure 1 and the 
average accretion disk My is shown by the thick gray line (Nelemans et al. 
2004). The system age (time since coming into contact) is indicated in the 
topmost panel by small circles spaced evenly in log(t/yr) and demonstrates 
that a cool, more compact donor leads to an older system at a given P^f,. In 
the Teff plot we delineate the DB instability strip (where isolated heUum rich 
WDs pulsate; Beauchamp et al. 1999). Constraints on T^ff for observed sys- 
tems are shown by either lower Umits or circles. (See Table 1 for a summary 
of measurements.) Some of the systems in the 20-40 minute range probably 
have a disk contribution even in quiescence, giving a low T^ff different from 
that of the WD surface. GP Com and CE-315, however, have both quies- 
cent disks and measured distances, giving T^fi and My directly from the WD. 
Comparing these to our models we find that both of these systems appear to 
have hot donors. 



percentile) gives a good maximum Af, which for forb > 40 
min can be fit by Af « 5 X 10Vr(^'orb/40min)4 3\ and thus 
yields the relation for the coldest possible O.65M0 WD at 
these orbital periods, L„„„ « 1 .5 x 1O~^L0(4O min/Z'orb)'*'^- 

3. OBSERVABILITY OF THE ACCRETING WHITE DWARF 

In order to assess the likelihood of observing the accreting 
WD, we must compare its luminosity to that from active ac- 
cretion. When the WD is near equilibrium or simply emitting 
the compressional heating luminosity, L ~ kTc (M) / ntp wiU be 
much less than the accretion luminosity, Lace ~ GM{M) /R, 
since kTc ^ GMnip/R. Seeing the hot WD then depends on 
either geometry (for example, an edge-on systmem like the 
newly discovered eclipsing binary SDSS J0926H-3624; Ander- 
son et al. 2005) or probes of different spectral regions. Four 
classes among observed systems are delineated in Table 1 by 
horizontal lines. The shortest period systems are in a phase 
where the accretion stream directly hits the WD (e.g. Marsh 
& Steeghs 2002; Israel et al. 2002), providing the first oppor- 
tunity. In this case there is no disc and there are times when 
the direct impact location is out of sight, providing an excel- 
lent chance to see the (very hot) accreting WD. At slightly 
longer Poib, around 20 minutes, all systems are expected to 



have a consistently bright accretion disc in the optical/UV that 
would energetically dominate the WD (e.g. Nelemans et al. 
2004) as noted above. Indeed, Nelemans et al. (2004) showed 
that a simple accretion disc model at 15 minutes would give 
My « 7 and temperatures of « 30,000 K, swamping the hot 
WD. This model is shown by the thick grey line on the bottom 
panel of Figure 2. 

The third group of AM CVn systems, with again slightly 
longer orbital periods, show large brightness variations due to 
a thermal instability in the disc (Tsugawa & Osaki 1997). Al- 
though at this (M) the disc absolute magnitude is expected to 
be brighter than the WD, the low state is a time when the ac- 
creting WD can be detected, just as in the hydrogen accreting 
CVs (Sion 1991; Townsley & Bildsten 2003). 

The last group, with the longest orbital periods, again show 
hardly any brightness variations. In the context of the thermal 
instability interpretation, these systems could be on the lower, 
stable, branch of the instabihty curve and thus would be cool 
discs, consisting mainly of neutral helium. The steadily ac- 
creting luminosity is only '--^10 times larger than the cooling 
WD (see §2), enhancing the prospects of seeing the accret- 
ing WD relative to shorter orbital periods, where the contrast 
is greater. In addition, the disc temperature will generally be 
cooler than the WD as shown in the bottom panel of Figure 2. 

3.1. Non-Radial Oscillations of the Accreting WD 

In addition to directly seeing the re-heated WD in the AM 
CVn binaries, there are also opportunities for seismologi- 
cal tests as the WD passes through the non-radial pulsation 
instabihty strip (Teg = 22,400-27,800 K for isolated DBs; 
Beauchamp et al. 1999), which we plot in Figure 2. The pul- 
sation periods are 200-1000 seconds (see Handler 2003 for a 
recent summary) and are typically / = 1 , 2 g-modes of high 
radial order (e.g. Kotak et al. 2003). 

We do not know, a priori, that the instability strip wiU be the 
same for a He accreting WD. The clear structural difference in 
the accreting WD is the mass of the helium envelope, which 
will likely exceed O.O5M0 and be seismically measurable. It 
seems unlikely that the thicker He envelope would alter the 
instabihty strip, since the relevant "pump" for the mode exci- 
tation occurs at low pressures (see Gautschy & Althaus 2002 
for a recent study). However, active accretion onto the outer 
atmosphere might make a difference, especially if the convec- 
tion zone (responsible for the mode excitation; e.g. Brickhill 
1991) is altered. It also seems possible that the mode pump 
could become active during the cooUng of the outer WD lay- 
ers (Piro, Arras & Bildsten 2005) following an accretion disk 
outburst. Indeed, Wood et al. (1987) reported excess power 
at periods of 100-200 seconds when CR Boo had an optical 
magnitude of V = 15.4, which might well be the passing of 
the accreting WD through the instabihty strip following an 
accretion outburst. 

The most dramatic alteration of the g-mode frequencies 
would be from WD rotation at periods comparable to the 
g-mode periods, as the WD has accumulated a significant 
amount of mass and angular momentum since the last mass 
ejecting He novae event. This angular momentum either re- 
sides solely in the He layer or has been shared with the under- 
lying C/0 WD in the ^ 10** years spent at the Porb « 20-40 
minutes where Jeff is in the isolated DB instability strip (see 
Figure 2). Motivated by the likelihood that very weak internal 
magnetic fields can easily mediate angular momentum trans- 
port on this timescale (e.g. Spruit 2002), we presume that all 
accreted angular momentum after the last He novae event is 
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TABLE 1 

AM CVn System Parameters and Inferred Accreting WD Properties 



System 


Porb 






Teff 


V Range 


Mv 


Ref 


Name 


(min) 


(M0 yr-l) 




(1000 K) 








RX J0806 


5.35 


1.7 X IQ-^ 


0.125 


> 40 


21.1 




1 


V407 Vul 


9.49 


7.0 X 10"' 


0.068 


>20 


19.9 




2 


ES Cat 


10.3 


4.4 X W 


0.062 


130±10 


16.9 




3,17 


AM CVn 


17.1 


3.7 X 10^"^ 


0.035 




14.0 


4.8 


4 


HP Lib 


18.4 


3.8 X 10-10 


0.032 




13.6 


7.0 


4 


CRBoo 


24.5 


7.6 X 10-11 


0.023 


12 


13.5-17.5 


5.7-9.7 


4,5,6, 7 


KL Dra 


25.0 


4.3 X 10-11 


0.022 




16.8-20 




8 


V803 Cen 


26.9 


7.6 X 10-11 


0.020 


14 


13.2-17.4 


5.1-9.3 


9 


CPEri 


28.4 


3.5 X 10-" 


0.019 


17 


16.5-19.7 




10,19 


SDSSJ0926+3624 


28.3 


3.5 X 10-11 


0.019 


18 


18-19 




20 


2003aw 


33.8 


1.4 X 10-11 


0.015 


17 


16.5-20.3 




11,18 


SDSSJ1240-01 


37.4 


5.6 X IQ-i^ 


0.013 


17 


19.6 




12 


GPCom 


46.5 


2.2 X 10-12 


0.010 


11 


15.7-16 


11.3-11.7 


4,13,14 


CE-315 


65.1 


3.4 X 10"" 


0.006 


9.5 


17.6 


13.2 


15,16 



References. — 1 Israel et al., 2002; 2 Ramsay et al., 2002; 3 Warner & Woudt, 2002; 4 Groot et al., 2005, in prep.; 5 Wood et al., 1987; 6 Provencal et 
al., 1997; 7 Patterson et al., 1997; 8 Wood et al., 2002; 9 Patterson et al., 2000; 10 Abbott et al., 1992; 11 Warner & Woudt, 2003; 12 Roelofs et al., 2005a; 13 
Thorstensen, 2003; 14 Warner, 1995; 15 Thorstensen, Priv. Comm.; 16 Ruiz et al. 2001; 17 Espaillat et al. 2005; 18 Roelofs et al, 2005b; 19 Sion, Priv. Comm.; 
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shared with the whole WD, yielding a lower limit to the ro- 
tation frequency of the He shell, where much of the g-mode 
energy resides. 

We assume that helium accretes onto the WD with the 
specific angular momentum of the inner edge of the accre- 
tion disk, {GMRy/^, and that the amount of accreted matter, 
AM <C M. The final WD spin frequency is then 27r/Pspin = 
(GMRy/^AM/I, where / ^ Q.IMR^ is the WD moment of 
inertia, yielding 27r/Pspin = 5AM(G/MR^y/^, 

f 0.05MQ \f M \''^ ( R 
npin « 50 s j \Q.mQ J V9 X lO^cm ) ' 

(2) 

which reaches about 100 seconds when M = \ .2Mq. Hence, 
it appears very likely that the WD surface layer (and maybe 
even core) is rotating rapidly enough (i.e. Pspin < ^mode) that 
the seismology will be dramatically altered from the isolated 
DB case. The discovery and study of a DB pulsator in an AM 
CVn binary would provide an unprecendented opportunity to 
measure the accumulated He layer mass on the C/O WD and 
the angular momentum transport mechanisms within it. The 
most probable orbital periods for this range from 20 to 35 
minutes. 

The rapid surface rotation will broaden any spectral fea- 
tures originating there, as the resulting minimum equatorial 
rotational velocity is w 1250 km s"'(AM/O.O5M0), nearly 
independent of M. The narrow, central line emission (i.e. 
spikes) seen in He lines of GP Com, CE 315 and SDSS1240 
(Marsh 1999, Ruiz et al, 2001, Roelofs et al, 2005a) are gen- 
erally believed to originate on the accreting WDs. However, 
their widths are much less than this estimate, suggesting either 
a different origin or an effective loss of angular momentum 
from the rapidly spinning accretor. 

3.2. Atmospheric Properties 

Another difference with single DB WDs is that is 
present in the accreting material at a mass fraction of 1-2% 
(due to CNO burning in the progenitor star of the He donor). 
This ^''N would rapidly sink out of the WD photosphere were 



it not for the constant mixing provided by the underlying 
convection zone of mass Mconv It is thus the sedimentation 
time, fsed, at the base of the convection zone that sets the 
relevant timescale. Paquette et al. (1986) calculated 4ed for 
'"'N in pure He atmospheres of isolated WDs, and defined 
an Msed = A#conv/fsed, bclow which the '"^N in the convective 
envelope will be less abundant than in the accreting mate- 
rial. For M = 0.6Mq WDs, Msed ~ IQ-'^Meyr"'. Though 
^conv decreases rapidly with increasing M (see Fontaine & 
Van Horn 1976), so does the sedimentation time, making the 
actual value of Msed not too dependent on M for the range of 
Tes relevant here (see Figure 14 in Paquette et al. 1986). More 
recent calculations with updated opacilies (MacDonald, Her- 
nanz, & Jose 1998) found larger Mconv values by factors of 
4-6 from Paquette et al. (1986). Hence, for PqA < 40 minutes, 
the atmospheric N abundance will be comparable to that in 
the accreted material. At longer orbital periods, depending on 
the state of the donor, it seems plausible that (M) will be less 
than Msed> and the '"^N will tend to sediment. Detection and 
study of atmospheric N (i.e. not from the accretion disk) in 
GP Com, CE-315 and SDSS J1240 could thus constrain the 
sedimentation physics. 

The abundant presence of '''N may simplify the spectral 
modeling of these WDs relative to what has been done for 
pure He atmospheres (see Bergeron, Saumon & Wesemael 
1995a). This is because the low first ionization potential of 
^"^N will provide free electrons more readily than the He. This 
may well allow the He" free-free opacity to more readily dom- 
inate for the relevant Teff's, but only calculations will tell. This 
could also alter the location of the pulsational instability strip. 

4. COMPARISON TO OBSERVATIONS 

We tabulate relevant observed and inferred properties of the 
AM CVn binaries in Table 1. The minimum donor mass, 
Mc,min, for cach system is the mass of the Deloye et al. (2005) 
fully degenerate donor model that fills the Roche lobe at 
the system's /orb- The minimum (M)min at fixed /orb oc- 
curs when Mc = Mc,min (Deloye et al 2005). We thus take 
A^cmin and calculate the tabulated values of (M)min as fol- 
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lows. For systems in which the mass ratio has been inferred 
(all the AM CVn binaries except RX J0806, V407 Vul, ES 
Cet, and SDSS J0926+3624), we use it and Me^min to find M 
and thus (M)mi„. When there is no mass ratio constraint, we 
assume that (M)min is given by the minimum possible value, 
which occurs at a mass ratio of ~ 0.3 (this determined the 
entries for RX J0806, V407 Vul, and ES Cet), except for 
SDSS J0926+3624, for which we use the mass ratio of CP 
Eri. 

Our discussion in the previous sections outlined the oppor- 
tunities for detecting the hot WD and we now discuss in more 
detail the implications of our calculations in regards to di- 
rectly observing the WD in these binaries. 

4. 1 . Short-period systems 

The evolutionary state of the two shortest period systems, 
RX J0806.3+1527 (Porb = 5.36 min) and RX J1914.4+2456 
(=V407 Vul, Porb = 9.48 min) is uncertain and different mod- 
els have been proposed (see e.g. Wu et al. 2002; Ramsay et 
al. 2002; Marsh & Steeghs, 2002; Marsh & Nelemans, 2005; 
DairOsso et al. 2005; Willems & Kalogera 2005). Marsh & 
Steeghs (2002) suggested that V407 Vul is a direct-impact ac- 
cretor where the accretion stream directly strikes the accreting 
WD. Noting the earlier work of Townsley & Bildsten (2002), 
they suggested that much of the optical light could originate 
from the hot, accreting WD. 

Our high (M) calculations allow us to address the direct im- 
pact scenario. Even though the accreting material arrives over 
a limited part of the WD, we expect complete spreading by the 
time it has reached depths relevant for compressional heating 
(Piro & Bildsten 2004), thus allowing our calculation to apply. 
At the minimum (M) » lO^^M© yr"^ expected for V407 Vul 
(see Table 1), we predict that L w L©. The predicted range of 
Teff w 50,000- 100,000 K and absolute magnitude. My « 8. 
The analysis of V407 Vul is complicated by the fact that the 
optical light is dominated by a G9V star (Steeghs et al., in 
prep). However, the variable part of the light clearly shows 
a blue color, with implied effective temperature > 20,000 K 
(see Table 1), consistent with our calculations. The estimated 
dereddened visual magnitude is 16.9 (Ramsay et al., 2002) 
which includes a substantial G-star contribution. Combined 
with the estimated My k, 8, this suggests a distance greater 
than 500 pc if V407 Vul is a direct impact accretor. 

For RX J0806.3-I-1527 (Porb = 5.36 min), the situation is 
quite different. The shorter orbital period implies a luminos- 
ity of « 44L0, Teff « 140,000K and My ~ 4.7 for the accre- 
tor. The estimated Tf.fi agrees with the observed very blue 
spectrum, for which Israel et al. (2002) find a lower limit 
of 40,000 K. The apparent magnitude of 21.1 would put the 
system at 20 kpc, 8 kpc above the Galactic plane. A similar 
distance is needed to interpret the rather low X-ray flux (Israel 
et al. 2002) as being generated by (M). 

It was recently suggested that the low X-ray flux could be 
reconciled with much smaller distances due to a temporar- 
ily lower M <C (M) or due to most of the accretion lumi- 
nosity being emitted in the UV (Marsh & Nelemans 2005; 
Willems & Kalogera, 2005). However, our calculations show 
that not only the X-ray luminosity, but also the optical flux 
tracks (M), so that the latter suggestion would still produce a 
much too bright WD. The compressional heating luminosity 
would track a drop in M on the thermal time of the accreted 
layer, r ~ cpT„(0.05Mq)/L ~ lO''yr((M)/lO-^M0 yr"')-" '' 
where is the maximum temperature, and L is given by 



equation (1). This timescale is much longer than the timescale 
for which a temporarily lower M can be sustained (~ 100 yrs. 
Marsh & Nelemans 2005), leading to a bright WD even at 
low M. Therefore, for RX J0806.3H-1527 to be a direct impact 
accretor, the only solution seems to be a very large distance. 

The third system in the table, ES Cet, is clearly a helium ac- 
creting WD. Espaillat et al. (2005) suggest it is a direct impact 
accretor and find Jeff « 130,000 K and radius /? » 4.5 x 10^ 
cm for a distance of 300 pc. They interpret this as the size of 
the impact spot, but our calculations suggest that the whole 
WD should be very hot. This observed temperature is in 
the expected range, and the radius either suggests a massive 
(I.IM0) WD or a less massive WD at a larger distance than 
they assumed. 

4.2. Outbursting sources 

For the outbursting systems the spectra are dominated by 
the accretion disc when they are in their high state. The ab- 
solute magnitudes of the disc are at least 6-8, much brighter 
than the expected absolute magnitudes of the accreting WDs 
of about 1 1 (Figure 2). However, in the low state, the WD 
can contribute significantly to the spectra, potentially produc- 
ing the observed continuum. Between 20 and 30 min, the ex- 
pected disc absolute magnitude drops, while the WD My stays 
roughly constant. Interestingly, the magnitude range observed 
in these outbursting systems decreases with increasing orbital 
period. While CR Boo has a range of 5 magnitudes, V803 Cen 
is less, at 4.2, and CP Eri only changes by 3.2 magnitudes. 

To check if the low- state spectra could be dominated by the 
WD, we estimate the temperature from the shape of the 4000- 
6000 A continuum from the few spectra that are available. We 
assume a black-body spectrum, which reasonably produces 
the broad band colors of theoretical DB white dwarf spectra 
(Bergeron, Wesemael & Beauchamp 1995b). Our tempera- 
ture determinations in Table 1 are accurate at about the 10% 
leveP. For CR Boo and V803 Cen, the absolute magnitudes 
are also known (Groot et al. in prep, see Table 1). 

For CR Boo and V803 Cen, the temperatures are below the 
expected WD temperatures and indeed the absolute magni- 
tudes in the low state are about 1 magnitude above the pre- 
dicted WD brightness (see Figure 2), suggesting a significant 
contribution from the disc (which is cooler than the WD, but 
of larger area), even in the low-state. Using the slope of spec- 
trum in the low state, we found a temperature of 30,000K 
for CP Eri, whereas a recent UV spectrum suggests 17,000K 
(Sion, priv. comm.). We have to be careful with these low- 
state estimates, though, as a second low-state spectrum of CR 
Boo (Patterson et al. 1997) shows a very different continuum 
shape, suggesting a cool or hybrid cool plus hot spectrum. 

The recently discovered eclipsing AM CVn binary, 
SDSS J0926H-3624, at Porb = 28.3 min (Anderson et al. 2005) 
provides an opportunity for probing the distinct contributions 
of the hot WD and the disk. Anderson et al. (2005) report 
that the « 1 minute duration of the ecUpse is consistent with 
the 0.02M© secondary transiting a hot spot or the inner disk. 
It is also consistent with an eclipse of the hot WD. Fast pho- 
tometry and phase resolved spectroscopy could confirm that 
the WD contributes and would tie the eclipse minimum to the 

3 For Teff < 17000 K the black body colors overestimate Teff by 500-1000 
K, while above that temperature they underestimate the temperature by a sim- 
ilar amount. In addition, in the region considered here, the change in the 
slope of the black body spectrum decreases with temperature above ~ 15000 
K, adding to the inaccuracy at higher temperatures 
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WD rather than the disk, yielding a more robust measurement 
of the orbital period change expected from gravitational wave 
driven mass transfer (Anderson et al. 2005). 

4.3. Long-period systems 

For GP Com and CE 315 we derived the temperatures in 
the same way as described above and interpret the temper- 
ature as that of the accreting WD. For GP Com, using the 
known distance (d = 67 pc; Thorstensen 2003), we can de- 
rive the emitting area. Thorstensen (2003) gives a typical ob- 
served V magnitude of 16. 1, but Warner (1995) lists the range 
15.7-16.0. We therefore use V=15. 7-16.1 as the observed 
range. Together with m-M = 4.2 ± 0.2 (Thorstensen 2003) 
this gives My = 1 1 .7 ± 0.4. From Bergeron et al. (1995b) we 
use the bolometric correction for a 1 1,000 K DB model, being 
-0.57. This gives a bolometric magnitude Mboi = 1 1 . 13 ± 0.4, 
i.e. a luminosity of 2.78 x IQT^Lq with an error of 44%. To- 
gether with the effective temperature, this gives the size of the 
emitting region of 10'^ cm, i.e. a typical WD radius! 

Using the same argument for CE 315, with a distance of 77 
pc (Thorstensen, priv. comm.), using a effective temperature 
of 9500 K and a bolometric correction of -0.34, we find an 
emitting region of 6 x 10** cm, implying a more massive WD 
of « O.9M0. Indeed, its absolute magnitude falls close to 
the track for a 1.0 Mq accretor (Fig. 2). Interestingly, two 
systems have been recently discovered that, in their low state, 
clearly show the broad absorption lines of a DB white dwarf, 
with estimated temperature of 17,000 K (SDSS J1240-01 and 
SN2003aw, Roelofs et al. 2005a; Warner & Woudt, 2003, 
Roelofs et al, 2005b). Based on an absolute magnitude of 
1 1 (Fig. 2), these systems should have distances of 735 pc 
(SN2003aw with V = 20.3) and 525 pc (SDSS J1240-01, with 
V= 19.6). 

5. SUMMARY AND FUTURE WORK 

We have shown (both from theory and observation) that 
thermal emission from the accreting WD within an AM CVn 
binary is an important contributor to the observed optical 
emission. For the direct impact systems with the shortest or- 
bital periods, the hot WD is visible when the hot spot is out 
of the line of sight. At longer orbital periods, the WD can ei- 
ther be seen when the disk is in a low state, or at even longer 
orbital periods (PqA > 40 min), the WD always dominates the 
accretion luminosity in the optical bandpass. 

We compared our theoretical results to the temperature and 
brightness of observed systems. The accreting WD is the 
dominant contributor to the optical light for the direct im- 
pact accretors and constrains their distances. In particular, 
the expected brightness of the WD is so large for RX J0806 



that it can only be in the direct impact phase if it is located 
at a distance of 20 kpc, far outside the Galactic disc. For 
Porb = 20 - 40min, our predicted absolute magnitude of the 
WD is My « 1 1 . This allows for estimates of the distance 
to the two newly discovered AM CVn systems (2003aw and 
SDSSJ1240-01) in which the DB white dwarf is visible. The 
two outbursting systems (CR Boo and V803 Cen) with paral- 
laxes have My « 10 in quiescence, 1 magnitude brighter than 
our estimate for the WD. This suggests that even in the low 
state there is a disc contribution. For GP Com and CE-315, 
a combination of estimated effective temperature and known 
distance suggests that the optical continuum emission is dom- 
inated by emission from the WD, in agreement with our cal- 
culations. For these two systems we find consistency with a 
hot (high-entropy) donor. Whether this donor was born hot 
and has evolved at fairly constant entropy or is heated by on- 
going processes is unknown but could be addressed by future 
theoretical work. 

Our successful comparisons to observations gives us con- 
fidence to apply our work to predict just how faint these bi- 
naries can become. For example, either mass WD reaches 
My = 15 at 5 Gyrs, and at 10 Gyrs, the absolute magnitudes 
are My > 20. This will allow for meaningful constraints 
on the population density of such binaries from deep optical 
searches, both in the field (e.g. the 3 additional candidates 
found by SDSS, Anderson et. al. 2005) and in Globular Clus- 
ters. 

Much work remains to be done on the theoretical side, the 
most important of which is to resolve better the nature and 
occurrence rate of He flashes at early stages in the evolution 
of these binaries. In addition to better constraining the re- 
maining He mass on the accreting WD at late times, it would 
also allow us to predict the He novae rate and potentially con- 
strain the population from the measurements of such events. 
Perhaps the most dramatic development would be the secure 
detection of non-radial oscillations in these WDs, allowing for 
measurements of both the WD rotation rate and accumulated 
He layer mass. 
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